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Introduction 


Tfiis report covers progress and activity on NASA Grant NS6-2016 
from !• April 1975 through 30 September 1975. Effort during this time 
period had three objectives: 

(a) Final development of a static aeroelastic analysis 
package for oblique wings in the subsonic flight 
regime. 

(b) Improvement of the subsonic flutter analysis package 
with three-dimensional unsteady aerodynamics. 

(c) Investigation of flutter analysis techniques suitable 
for the supersonic flight regimes. 

Results of these efforts are discussed below. 

The documentation and source deck for the static aeroelastic analysis 

computer program with the acronym AIRLOD was delivered to NASA Ames Research 

Center in May 1975. This program was forwarded to LTV - Dallas by NASA/Ames. 

Checkout of this program by LTV has proven satisfactory; a few minor errors 

were uncovered and corrected. Program extensions and improvements have 

been suggested, but these improvements are being held in abeyance because 

of more pressing problems. 

Subsonic Flutter Analysis 

Several improvements have been made in the subsonic flutter analysis 
package. Chief among these improvements is an improved, more accurate 
mass matrix formulation and capability to handle three rigid-body degrees 
of freedom. In addition, ar. expanded modal analysis capability was added 
to the program package. Checkout of the analysis package has progressed 
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to the.potot urtiere a hi(^ conftdep.ee level in toe progrem has been attotned^^^^^ 

A preliminary copy of this progran package has been delivered to NASA/Anes 
for use on the Ames CWI 7600. The program was originally developed for 
the I.B.M. 370/158 machine. 

The AIRLOD computer program and the subsonic flutter analysis computer 
programs have been used to generate a great deal of new numerical data. 

This data has formed the basis of two papers written during the summer 
and submitted to the AIAA Journal and AIAA Journal of Aircraft . Ctotes 
of these papers have been sent to the grant nranitor and ar^ attached to 
this report. 

Principal Theoretical Results 

Principal findings thus far have been that at least too types of flutter 
instability can occur in the subsonic flight regime. The first of these 
types of instability is associated v.ith a low value of the ratio of wing 
roll moment of inertia to fuselage roll moment of inertia. This mode of 
instability is characterized by low reduced frequencies (k) at flutter 
(k ® wb/V, where w = frequency, b = wing semi -chord and V is airspeed). 
This type instability involves primarily a coupling between wing elastic 
bending and fuselage rigid roll; this instability resembles a vehicle in- 
stability. 

The second type of flutter instability involves classical bending- 
torsion flutter with i slight amount of rigid roll coupling. This mode 
of instab' 11 ty occurs at a relatively high reduced frequency (of the order 
of k = 0.25) and is characterized by larger values of the ratio of wing 
roll moment of inertia to fuselage moment of inertia. In general, if the 


transition from bending- tors ion flutter to bwwiing'^rolV flutter can b* 
precluded, the flutter speed of the aircraft can be kept at yalues_far 
above the wing's clamped divergence speed. Research is cmtinuing to 
discover the mechanism which triggers this change. 

Although primary emphasis has been thus far plaod on ^e effect 
of roll on the flutter of oblique wings, research has been initiated on 
the effect of fuselage pitch inertia on stability. The objective of this 
research will be to define critical parameters which may adversely couple 
rigid body pitch to elastic deformation. 

Transient Response 

It should be noted that the aerodynamic influence coefficients generated 
in the flutter program also may be used to solve for the transient re- 
sponse of the vehicle at low reduced frequency. These low reduced fre- 
quencies correspond to a slowly maneuvering vehicle. Such a study of 
dynamic stability and control is not turrently within the scope of the present 
grant and is deferred to some future time. However, subsonic response of 
the aircraft at speeds below the critical speed is of great importance. 

Supersonic Flutter Analysis 

The supersonic flight regiiria with subsonic leading edges is, for 
flutter, next in importance after the high subsonic regime. Unfortunately, 
no one unified method seems available for the generation of supersonic 
aerodynamic influence coefficients. The particular difficulty which 
arises in the case o; ;he subsonic leading edge is that the region between 
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the foremost Mach cone and the leading edges (the so-called diaphragm 
region) complicates the solution for velocity potential formulations. 

The Mach box formulation of the supersonic flow problem presents 
one way of solving for the aerodynamic influence coefficients. This 
method has been found to have several deficiencies, the principal 
shortcomings being inaccuracy at low supersonic Mach numbers (M < 1.2) 
and the computational time necessary for some planform configurations. 

In addition, available computer codes are restricted to analyzing sym- 
metrical planforms. A literature search has uncovered three potentially 
worthwhile methods for solving the supersonic problem; these methods are 
discussec below. 

Available Solution Methods for Supersonic Flutter Analysis 

A computer program is available through COSMIC which can be used to 
predict flutter at supersonic speeds. This program is called LAR-10199 
and was developed by E. C. Yates at Langley Research Center. The 
analysis method used in the program is based upon so-called modified strip 
theory. A separate analysis of the static, rigid wing must be available 
to predict the spanwise distribution of lift curve slope and the chord- 
wise position of the aerodynamic center (AC). This latter quantity is 
extremely important to flutter analysis. It has been noted that the 
failure of linearized supersonic theory to accurately predict the AC 
position leads to overly unconservative flutter predictions in many cases. 
However, Yates shows that with the proper prediction of the static AC, 
good results may be achieved. 


A second method of supersonic flutter prediction Is the use of piston 
theory. While piston theory is commonly used to predict flutter when 
M » 1 or M k » 1, a correction has been suggested to give agreement with 
the second-order quasi-steady supersonic theory of Van Dyke. This cor- 
rection is thought to extend the validity of piston theory to lower 
supersonic Mach nunters and lower reduced frequency. 

The Mach box method is also available for use in the flutter analysis. 
It appears, however, that a reprogramming of the method is necessary to 
account for asymmetry of the wing. This roprogramning effort does not 
need to be as extensive as existing Mach box programs since existing 
programs usually concern themselves not only with the wing pressure dis- 
tribution but also with wing-tail interference effects. This latter 
effect is of no concern to the present analysis. An investigation of 
the time required for this programming effort is currently underway. 

In addition to the above methods, other techniques have been suggested 
for calculating supersonic AIC's. Prime among these are the acceleration 
potential kernel function techniques. A literature survey of related 
papers is currently being conducted. 

Planned Approach to Supersonic Flutter Analysis 

Because of the disparity > f analytical results often encountered in 
supersonic flutter analysis, the author believes it to be prudent to pursue 
several courses of study simultaneously. As a temporary measure, the 
program LAR-10199 should be obtained and made operational as quickly as 
possible. At the least, this method can identify, in a qualitative 



manner, the flutter behavior which may be encountered in the supersonic 
flight of the oblique wing. 

A listing of a computer program to conqiute piston theory AIC's with 
thickness effects and sweep corrections included is available at V.P.I. 

This program has the capability of generating aerodynamic data for a flutter 
program. An attempt will be made to make this program operational. 

Finally, the effort necessary to generate a simple Mach box routine 
will be assessed, together with an investigation of newer supersonic AIG 
generation methods. The rationale behind the study of these methods is 
to generate and compare data generated by different theories and to coiir- 
pare differences and similarities of results. 

Continuing Subsonic Flutter Work 

A substantial amount of subsonic flutter investigations remain to 
be done. These studies include the effect of pitch inertia on flutter 
and the effect of shifting the wing c.g. relative to the aerodynamic 
center and the elastic axis. These efforts will continue; a graduate 
student is presently doing graduate work in this area. 

Travel 

During the time period covered, a trip was made to the AlAA Structures/ 
Structural Dynamics Meeting held in May 1975 in Denver, Colorado. In 
addition a visit to NASA/ Ames Research Center was made on 25>26 September 
1975. During this visit the latest results of subsonic oblique wing 
flutter studies were presented. 



FLUTTER OF ASYMMETRICALLY SWEPT WINGS 


by 

Terrence A. Weisshaar* 
and 

J.B. Crittenden'*’ 

Virginia Polytechnic Intcltute 
and State University 
Blacksburg, Virginia 


This research was supported by NASA Ames Research Center, Moffett Field, 
California under NASA Grant NSG-2016. 

Index categories: Aircraft Handling, Stability and Control; Aeroelastlclty ; 

Aircraft Structural Design. 

^Assistant Professor, Aerospace and Ocean Engineering Department, Meadier 
AIAA. 

'^'Research Assistant, Aerospace and Ocean Engineering Department. 


FRBCEDtNO PAGE BLANK WOI 


ABSTRACT 


Recent Interest In asymmetrically swept, oi' oblique, wings has raised 
fundamental questions about the flutter characteristics of such wings. This 
paper presents two formulations of the oblique wing flutter problem; one forsaf 
latlon allows wing bending deformations and the rigid body roll degree of freedom 
while the second formulation Includes bending-torsional deformation and roll 
degrees of freedom. Flutter is found to occur in two basic inodes. The first 
mode is associated with bending-roll coupling and occurs at low reduced frequency 
values. The other instability mode is primarily one of classical bending-torsion 
with negligible roll coupling; this mode occurs at much higher reduced fre- 
quencies. The occurrence of bending-roll coupling mode leads to lower flutter 
speeds while the bending-torsion mode is associated with higher flutter speeds. 
The ratio of the wing mass moment of Inertia in roll to the fuselage moment of 
inertia evidently plays a major role in the determination of which of the two 


instabilities is critical. 
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Nowenc'lature 


[a] * flexibility matrix for clamped fuselage vring 
c » wing chord measured perpendicular to elastic axis (Fig. 1) 
c = wing chord measured parallel to the free stream direction 
Cl, = 2-dimensional sectional lift-curve slope 
g = structural damping parameter 

i = »cr 

I = wing roll moment of inertia at zero sweep 
0 

k = reduced frequency, kc/2Vjj or uc/2V 

V = airspeed 
V|r = flutter speed 

= airspeed normal to swept axis, = VcosA 
A = sweep angle 
p = air density 
b) = frequency of oscillation 


V 
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Introduction 


The recent interest in the use of an asyranetrically swept, high- 
aspect-ratio wing to achieve high lift-to-drag ratios has generated . 
interest in the aeroelastic stability characteristics of such a configu- 
ration. However, the undesirable static aeroelastic behaviw* of 
symmetrically swept forward wings has prompted some caution on the part 
of structural engineers towards the asymmetrical wing. As a result, con- 
siderable discussion of the merits of such a design and the potential 
weight penalties which might be incurred has occurred. Jones and Nisbet 
(Ref. 1) have presented data which tends to allay some misgivings about 
the aeroelastic stability ''f asymmetrically swept or oblique wings. Pro- 
minent among their findings is the discovery that the inclusion of the 
rigid body aircraft roll degree of freedom has a stabilizing effect on 
the aeroelastic stability of the wing, when compared to the stability of 
a similar, but clamped, wing. Their analytical results were obtained 
through the use of quasi-static aerodyanmic theory to represent the per- 
turbation lift forces generated by the harmonic motion of their idealized 
flexible model. 

This study seeks to explore in somewhat more detail than Ref. (1) 
the flutter behavior of asymmetrically swept or oblique wings; to accom- 
plish this task the results of two studies are presented. The first 
study examines the flutter behavior of an idealized, uniform-property 
wing in incompressible flow and at various asymmetrical angles to the 
flew. For this portion of the study, quasi-steady aerodynamic strip 
theory will be employed in the equations of motion; the Galerkin method 
will be used to solve these equations. 
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The second portion of the study entails the use of a more sophlsti* 
cated approach to the solution of the oblique wing flutter problem. This 
approach uses a finite-element* unsteady aerodynamic representation 
together with a multl-degree-of -freedom structural model to examine more 
closely and more accurately the flutter behavior of variable planform 
wings. 

From these studies, It will be shown that, at moderate sweep angles, 
the flutter speed of the wing may be lowered when compared with the flut- 
ter speed of the wing at zero sweep, addition, the shap' of the wing 
planform and the spanwise dist**1but1on of stiffness and weight will have 
a significant effect on the relation between flutter speed and sweep 
angle. 

Discussion 

The first part of this study Is concerned with the aeroelastic anal- 
ysis of a simplified oblique wing model, shown In Fig. 1. The impetus 
for such a study stems from the desirability of assessing the behavior of 
the flutter speed of the wing as it Is asymmetrically swept. This model 
repre5 ;nts a wing of uniform structural and aerodyanmic properties, asym- 
metrically swept at an angle jv to the flow. This high aspect ratio wing 
Is idealized as a beam with a straight elastic axis, free to roll about 
an axis poialle! to the flow. It is assumed that mass is distributed 
along this roll axis such that a mass moment of inertia, simulating 
the roll moment of inertia of the fuselage, appears concentrated there. 

To examine the aeroelastir stability of this model, assume that It 
is caused to undergo small ostilfatlons about a "wings level" static 
eqW librium portion. The stability of the subsequent motion can be 
determined by an examination of the character of this free vibration. 
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The structural oehavior of this wing can be modelled through the use of 
conventional Euler-Bernoulli beam theory. It is further assumed that the 
wing has no torsional flexibility so that only bending flexibility is 
important. The limits to the validity of this latter assumption will be 
discussed later in this paper. 

In Ref. 2 , Barmby and Cunningham discuss the flutter analysis of 
symmetrically swept wings through the use of aerodynamic strip theory and 
the Theordorsen functions. The present study neglects all the noncircu- 
latory aerodynamic terms in Ref. 2 but retains two of the circulatory 
terms. In addition, the free vibratory motion is assumed to take place 
at a value of reduced frequency k = o)C/2V^whi ch is so small that the flow 
ib quasi-steady. The circulatory aerodynamic terms retained are the term 
which corresponds to the familiar damping-in-roll and the term which 
arises from the angle of attack generated by bending deformations of a 
swept wing (Ref. 2). 

The assumptions about the behavior of this idealized irodel under- 
going small oscillations in the airstream lead to the following 
differential equation of motion for the elastic wing. 


.4. 


O I'l . Pf 3 ‘*1*1 . / _ __ rnc^A^ ^1^ 

m ^ + El ^ + (qcc^.cos A) 


tanA 


3t 


ay 


- mpycosA = 0 

where m = wing mass per unit length along the y-axis 

El = bending stiffness of cross-section perpendicular to y-axis 

q - freestream dynamic pressure 

' * 

\ W = wing upward deformation due to bending only 
t = time 

p = roll rate in radians per unit time 


( 1 ) 
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Nondimenslonalization of Eq. (1) yields the following equation. 




( 2 ) 


n = 0 

where {’) = differentiation with respect to time 
w = W/L 
n = y/L 

3 

X = qcc^^^L sinAcosA/EI 

The requirement that the sum of all roll moments generated by wing oscil- 
latory motion be equal to zero results in the additional equation: 


fif + IqCos^A^ P = ^gcc^^L^cos^A^ ^ J ^ ndn^ 
f 3 \ ■ 

+ I mL cosA 1 J wndn 

Z qcc, L^cos\ / dL\ 

■ 3 \V / 

qcCj^^L^ros^A . 

+ ^ J Wndn 

If we let (j) = p cosA then Eq. (3) may be written as 
^ + l(v^) ^ (vsinA) 

•w/.: 


ndn I tanA 


(3) 


/ Sw . . 3 /Iw 



(4) 


T/y-1 


1 . 

Wndn 


where Ij = + IqCOs A 

2 2 3 2 

I,^ = l-cos A = Y ml cos A 

wo 3 t 

2 3 * 

Y = qcc^^ L cos A/Iy 
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To solve Eqs. (2) and (4), the time dependency is eliminated by recoyni 
tion that the functions w(n,t) and 4{t) are separable such that 

w(n,t) = f(n)e*"^ 

Tt 


(5a) 


♦(t) " ♦€ 

Next, Eq. (2) is separated into two parts, one valid in the region -1 < n 
^ 0, the other valid in the region 0 < n ^1. Finally the resulting set 
of equations is solved approximately through use of Galerkin's method. A 
simple polynomial to use for such a solution is that shape obtained for 
uniform loading of a cantilever beam. In this case the function f(n) is 
approximated as: 


f(n) = 


B tf- 2 ii 3 , 4\ 

y (6n - 4n + n ) 


^ 1 


y (6n^ + 4n^ + n^) -1 - n - 0 


(6) 


r T 1 

a 


0 

|diil 1 

b 


0 

L’JJ| 

$ 
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where a and b are unknown constants, the Galerkin method leads to a set 
of three homogeneous algebraic equations, represented in matrix form as: 

(7) 

The coefficients d^^ are given in the Appendix to this paper. 

It is found that, in the absence of the roll freedom, the first 
natural frequency of vibration of the clamped winq, in vacuo , is pre- 
dicted by the Galerkin method to be 

0 )^ = 3.530 /117 ( 8 ) 

This compares with the exact solution (Ref, 3) 

0) = 3.518 [eT 


(9) 


Fo«;^the clamped wing, it is found that the sweptforward wing undergoes 
static divergence when r is zero. This occurs at a value of X equal to 
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6.40. The exact wolution gives a value of X for static divergence of 
6.33 (Ref. 4). 

If all the system parameters, such as El, A and V, are substituted 
into the expressions for d.., then the determinant of the matrix d.., 

I J I J 

written as A(d^j), can be used to find r through the relation: 

A(d^j) = 0 (10) 

With reference to Eqs. (5a, b), it is seen that if r is found to be a real 
number, then motion is aperiodic. A positive real value of r indicates 
aperiodic instability or static divergence. On the other hand, if r is 
found to be a complex number, motion is harmonic. If r = « + i« then the 
motion is periodic of frequency u. For negative values of «, the motion 
decays, but for positive values of « it grows with time. This latter 
situation corresponds to the dynamic aeroelastic instability commonly 
referred to as wing flutter. At the value r = iu, the system undergoes 
undamped oscillation and is said to be in neutral equilibrium. For a 
given set of system parameters, the airspeed V at which this occurs is 
called the flutter speed, Vp. 

The way in which the problem is presently formulated allows us to 
pick one system parameter as the unknown in Eq. (10). The maghitude of 
the complex number r is of no interest, but rather the value of velocity 
at which neutral stability occurs. For this reason, it is found to be 
advantageous to let r = iu in the expressions for d.. and to express these 

’ J 

coefficients in terms of and the parameter $, defined as 

e = x/x^.y (11a) 

V * 

where X^.^ = 32/5 = 6.40 (11b) 

The expressions for d.. in terms of these parameters may be found in the 

* vl 

Appendix. 


Given the systan physical parameters, the determinant in Eq. (10) 
may be expressed in terms of the independent variable 8. Collecting 
terms, the determinant is found to have a real part and an imaginary part 
given respectively by the expressions 



and 


+ fl - 6^ + 26 8R + iL 'I'dI “ 0 

L ^ 60 J 


(12b) 


+ ^ (1 + 8^/25) = 0 


where 


and 



(13a) 


(13b) 


The selection of a value of 8 which yields identical values of the 
ratio id/oJq in Eqs. (12a, b) completes the solution. From this value of 8, 
the flutter velocity may be obtained. 

The above solution procedure was implemented for a small model wing 
constructed of aluminum sheet with a constant thickness of 0.064 inches. 
The wing properties were taken to be: 

Material density = 0.101 Ib^^/in? 
c = 4 in. L = 20 in. Cj^^ * 2ir 

= 3 El = 874.0 Ib^-in. 

Using a sea level air density value, Eqs. (12a), (b) were solved numeri- 
cally using a Newton's method trial and error solution technique. The 

V * 

results of this analysis are shown in Fig. 2. 
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From Fig. it is seen that the flutter speed decreases as the wing 
is swept. For small values of A, the value of Vp greatly exceeds that of 
the clamped wing static divergence speed, Vq. However, as A increases, 
the critical speeds Vp and Vp draw closer together; at A = 90* they will 
coincide. As suggested by Jones and Nisbet in Ref. 1 , the moment of 
inertia ratio lyij plays a significant role in the flutter analysis of 
this asymmetrical wing. From the expression for it is clearly 

seen that this ratio tends to zero as A approaches 90**. It has been sug- 
gested that this mass moment of inertia ratio should be as large as 
possible to improve flutter performance. The results in Fig. 2 supi^ort 
this suggestion. 

Since one of the original assumptions of the present analysis was 
that the flutter instability occurs at relatively small values of reduced 
frequency k it is worthwhile to note the values of reduced frequency for 
which the instabilities in Fig. 2 occur. These numbers are listed in 
Table I. Although these reduced frequencies are reasonably small, the 
accuracy of these results is probably somewhat degraded by the quasi - 
steady flow assumption. 

The model just analysed is similar to, but not identical to, a 
series of models used by Papadales (Ref. 5) in wind tunnel experiments at 
Virginia Polytechnic Institute. Those experiments had as their primary 
objective the study of the static aeroelastic characteristics of clamped 
oblique wings. However, when those tests were completed, simple flutter 
tests were conducted on roll -free models. Although no attempt to take 
accurate data was made during these flutter demonstration tests, the 
ve^city magnitudes shown in ^fig. 2 correspond to the order of magnitude 
of the velocities observed in these demonstrations. In addition, for 



sweep angles greater than 15-20®, the primary mode of instability was 
observed to be a fundamental symmetrical bending mode coupled with a 
rigid body roll oscillation (Ref. 5). 

For sweep angles less than about 20®, the tests described in Ref. 5 
found a flutter mode which resembled a more conventional bending-torsion 
coupling with a slight degree of rigid body roll interaction. These 
observations, together with the desire to obtain a more accurate versa- 
tile analysis model, suggested the application of a more sophisticated 
analysis technique to the oblique wing flutter problem. It is to this 
analysis that we now turn our attention. 

Conventional flutter analysis of realistic aircraft employs assumed 
structural deflections or mode shapes together with generalized coordi- 
nates assigned to these mode shapes. An excellent discussion of modal 
and non-modal matrix methods of flutter analysis is given by Rodden in 
Ref. 6. In addition. Ref. 6 presents a succinct discussion of how to 
include free-free boundary conditions into the conventional restrained or 
clamped model. This latter discussion follows the development given in 
Ref. 7, but is more general. The highlights of Ref. 6 are reviewed here. 

To analyze the flutter behavior of a planform such as that shown in 
Fig. 4, it is necessary that the following items be taken into account: 

s. 

the distributed mass of the wing along the span; the variable bending and 
torsional stiffness along the span; and the unsteady, three-dimensional 
aerodynamic forces and moments associated with deformations caused by 
wing oscillations. With the assumption of simple harmonic motion at fre- 
quency bi, the classical matrix equation for flutter analysis, before the 
inclusion of assumed modes, ^s given by (Ref. 6): 
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In this equation the static flexibility matrix [a] has been divided by 
the factor (1 + ig) to account for the structural damping necessary to 
sustain simple harmonic motion. The elements of the vector {h} are 
actual elastic deflections and rotations at concrol points on the wing. 
The mass matrix [M] and aerodynamic influence coefficient matrix [C^] are 
both multiplied by the frequency squared. The coefficients of refer 
to the air density p, the reference semi -chord b^ and the wing semi -span 
s. The elements of [C^] are complex numbers and functions of Mach number 
and the local control point reduced frequency, k = ub/V, where b is the 
local semi -chord. With the formulation in Eq. 14, the unsteady aerodyna- 
mic forces enter into the problem, mathematically, as complex masses. 

The idealization of the wing structure as an assemblage of beams, 
each with a straight elastic axis, permits the use of conventional finite 
element structural analysis methods to describe the wing stiffness and 
flexibility. The reader is referred to Refs. 8 and 9 for discussions of 
this method. Similarly, the mass matrix may be formulated from finite 
element methods. The mass matrix must account for the fact that the wing 
shear centers may be offset from the wing chordwise location of the cen- 
ters of mass. Finally, to model the three-dimensional aerodynamic forces 
and moments, the doublet-lattice method (Ref. 10) was used. To employ 
this method, an existing computer program (Ref. 11) was employed to gen- 
erate aerodynamic influence coefficients. 

A computer program was written to calculate the matrices in Eq. 14 
using these methods. The free vibration modes for the clamped system 
are then used to reduce the size of the matrix equations. The free-free 
b^ndary conditions are then i introduced to "free" the clamped system 
described in Eq. 14; this allows rigid body roll Once these 

matrices have been formed, the eigenvalues and eigenvectors may be 
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found. Since the aerodynamic influence coefficients are functions of 
reduced frequency k and Mach number (in these studies, Mach number is 
zero), a set of eigenvalues and eigenvectors corresponding to each value 
of k is generated. The familiar V-g method (Ref. 3, pp. 565-568) is then 
used to find the value of velocity at which neutral stability occurs. 

To assess the effect of torsion and unsteady aerodynamics on the 
flutter analysis of the oblique wing, the uniform property aluminum wing 
was reanalyzed. The wing is considered to have the same structural pro- 
perties as before, but, in the present example, 6J is taken to be equal 
to 1346 lb-in. It should be noted that the flat, sheet-aluminum wing has 
a ratio of first bending to first torsion which is slightly higher than 
that common to conventional aircraft. 

The analysis of the constant property wing, including roll freedom 
and torsional flexibility and employing the doublet-lattice method was 
conducted with a sixty degree-of-freedom model. These sixty degrees of 
freedom were obtained by considering ten control points on each wing; 
each control point has pitch, plunge and bending rotation elastic degrees 
of freedom. This model was subsequently reduced to a twenty degree-of- 
freedom model by using the first twenty natural modes of the system. 

The results of this flutter analysis are displayed in Fig. 3 as 
ratios of the instability velocity (either flutter or divergence) to the 
velocity at which wing torsional divergence occurs at zero sweep; this 
latter velocity is denoted as V^q. 

In Fig. 3 the behavior of the wing when the fuselage is clamped is 
shown as the curve labelled Vq/Vjjq. With the fuselage clamped, instabi- 
lity is found to occur at a neduced frequency k = 0; this is divergence. 
\ * 

When roll freedom is allowed, and when I^/If = 3, a dynamic instability 
appears; this is flutter and is shown as the curve Vp/V^Q. Unlike the 



strip theory results, the flutter speed does not tend to Infinity as A 
tends to zero. In fact, the flutter speed reaches a maximum near 15* of 
sweep and then declines, approaching the clamped divergence speed at a 
sweep angle of 60*. From an examination of the mode of Instability in 
the analysis and from test observation, it is found that, from A * 0* 
until near A * 30*, the instability resembles bending-torsion flutter 
such as might be observed on symmetrically swept back wings; an increas- 
ing anrount of rigid body roll appears as A increases. 

As a further illustration of the flutter behavior of oblique wings, 
a non-uniform wing planform, constructed in the same manner as the uni- 
form property wing, was analyzed. This wing (Fig. 4) has a modified 
elliptical planform and is constructed to give, theoretically, the mini- 
mum induced drag for a given lift, given span and given root bending 
moment (Ref. 12). In this case, the wing-fuselage combination has a roll 
moment of inertia ratio Iq/I^: = 11.69. Fig. 5 shows tha stability behav- 
ior of the clamped and roll-free wings. While the decrease in divergence 
speed with increasing A for the clamped wing shown in Fig. 4 resembles 
that of the uniform property wing, the behavior of the flutter speed for 
the nonuniform wing is much different. Orice again, for large sweep 
angles, the decrease of flutter speed with sweep angle is seen in Fig. 5; 
however, the roll-free flutter speeds and clamped divergence speeds are 
more widely separated in Fig. 5 than in Fig. 3. 

To assess the importance of the roll moment of inertia, the aero- 
elastic stability of the nonuniform property wing shown in Fig. 4 is 
again studied. However, the roll moment of inertia of the fuselage is 
now increased by a factor of 'two. The results of this study are pre^r 
sented in Fig. 6 and are compared to those previously obtained using the 


smaller fuselage roll moment of inertia. Once again, the results are 
displayed as ratios of flutter speed to clamped divergence speed at zero 
sweep angle. 

The effect of increasing the fuselage roll moment of inertia is 
clearly seen in Fig. 6. The flutter speeds for both moment of inertia 
ratios are seen to be identical until about a 15 degree sweep angle. 

Near this point the flutter mode for the Iq/I^ = 5.85 wing becomes pre- 
dominantly one of fundamental bending with rigid body roll coupling. 

This is seen to depress the flutter speed as A increases. 

As a final example, consider the uniform property aluminum wing. 

This wing has been previously analyzed using quasi -steady strip theory 
aerodynamics and elastic bending degrees of freedom with roll coupling. 

It has also been analyzed with a bending-torsion model which used the 
doublet-lattice aerodynamic loads. For the present example, the value 
of the torsional stiffness GJ is chosen to be 10 times that of the exam- 
ple whose results were presented in Fig. 3. The results of the present 
study are shown in Fig. 7, where they are compared with those presented 
in Fig. 3. In Fig. 7, the designation "Uing 2" refers to a uniform 
property wing with properties identical to those of "Wing 1" except that 
wing sectional torsional stiffness of Wing 2 is ten times that of Wing 1. 

In Fig 7, the relation between flutter speed and A for Wing 2 is 
seen to have a discontinuity near A = 15*. Increasing the value of GJ is 
found to have a pronounced effect on flutter speed at moderate sweep 
angles, but has little effect on flutter at high sweep angles. The dis- 
continuity is believed to be caused by a sudden change in the' flutteir' 
mode near 15 degrees sweep! j For sweep angles beyond 15®, the results 


obtained for this wing with high torsional stiffpjss resemble those 
obtained with the bending model and strip theory airloads. The reduced 
frequencies at the onset of flutter of Wing 2 are displayed in Table II 
for several sweep angles. From this table, it is seen that the flutter, 
which occurs primarily as a roll-bending or flapping instability, occurs 
at relatively low reduced frequencier ;hert compared with the reduced fre- 
quencies which arise at the onset of bending-torsion flutter. Also, a 
comparison of the reduced frequencies in Tables I and II shows that the 
reduced frequencies at flutter in the two studies are comparable in mag- 
nitude. 

Conclusions 

Before summarizing the results of this paper and listing conclusions, 
certain features of the idealized models studied should be reviewed. 

These models were chosen for analysis because of past experience with 
wind tunnel tests. A constant thickness, sheet metal wing has a bending 
stiffness which is proportional to the wing chord measured perpendicular 
to the wing elastic axis; the torsional stiffness varies in a similar 
manner. This proportionality of the stiffness to the wing chord leads to 
bending and torsional stiffness distributions which are concave downward 
when plotted versus the spanwise coordinate. In actuality, the bending 
stiffness distribution which results from considerations of wing strength 
usually appears to have a concave upward distribution (cf. Ref. 3, p. 45). 

The wings considered in this study had wing sectional centers of 
mass coincident with the shear centers; thus, there was no elastic axis- 

c.g. offset. Dynamic coupling was either nonexistent, as in the case of 

' \ 

thVunii property wing, of minimal, as in the case of the variable 
property wing. This latter wing has a line of shear centers which is 
curved slightly forward when the wing is in its unswept position. 
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The combinations and permutations of the various parameters which 
affect the aeroelastic stability of an aircraft are seemingly endless. 
However, several conclusions may be drawn from the present studies at 
zero Mach number. Prominent among these conclusions Is that the Inclu- 
sion of the rigid body roll degree of freedom Into the flutter model 
causes the mode of Instability to change from an aperiodic Instability 
(divergence) to an oscillatory Instability (flutter). The degree to 
which the stability boundary Is modified depends upon the sweep angle A 
and the ratio of the moments of inertia in roll of the aircraft fuselage 
and the wing in its unswept position. 

If the wing instability appears as a coupling between the wing's 
fundamental bending mode and rigid body roll (the "flapping" mode), flut- 
ter speed is reduced as A increases. However, If the system p»ra~icters 
are such that flutter appears primarily as a bending-torsion coupling, 
the flutter spaed may actually increase as the wing is swept. If the 
wing can be either elastically or dynamically tailored, it may be possi- 
ble to avoid the flapping mode type of instability altogether. 

Topics V **ranting further investigation include: the effect of Mach 
number on oblique wing flutter; the significance of elastic axis - c.g. 
offset; and the effect of elastic tailorin" uT the wing. It is antici- 
pated that these and other studies will provide further insight into this 
unique aerodynamic design. 
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Appendix 

The application of Galerkin's method to Eqs. (2,3) results in the 


determinant in Eq. 10. The elements of the matrix |d^ 

jj are given below. 

<’ll ' + i«D + 1 + e 

(Al) 

d,2 = 0 

(A2) 

<•13 - 1 [fe)' - •“»] 

(A3) 

dj, =0 

(A4) 

<*22 ° • (l^) + *»0 + 1 - 6 

(A5) 

'*23 - - 1 [(^) - •“"] 

(A6) 

dj^ = u»^(13R/30) - i«<i(13/45) - (YsinA)(3/5) 

(A7) 

d 32 = - uj^(13R/30) + i(o<|»(13/45) - (YSinA)(3/5) 

(A8) 

*^33 ~ ~ ^ ^ io)ij»{2/3) 

(A9) 


The following definitions of terms are used in the above equations. 

D = (13/162)(xL/VsinA) t|» = yL/V 

R = (I^/Iy)cos^A X = qcC|^^L^sinAcosA/EI 

^ o 'n 

^ 3 = ^^^cr ~ * y ~ qcC|^^L cos'^A/Ij 
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Abstract 


A static aeroelastic phenomenon unique to an aircraft with asymnetHcally 
swept wings is disc ssed and the potential magnitude of its importance assessed. 
For this assessment, a simple formula is derived from the analysis of a highly 
idealized model. The validity of this formula is examined through the use of 
a more sophisticated numerical analysis. Among the results of this analysis 
are the following: aileron settings of a few degrees are sufficient to trim 

such aircraft in roll for Ig flight; the use of so-called built-in twist in 
the form of an initial wing anhedral provides an efficient alternative to 
aileron trim; if the wing is elastically tailored in a proper Tashion, it 
may be possible to design a wing whose elastic deformation under airloads 
provides a form of self- trim in roll at the cruise q of the aircraft, in 
which case no aileron input or anhedral is necessary. 


Introduction 


I . ., 1 .-.^ , 


The oblique wirg aircraft concept Is currently being explored as a possi* 
ble method of achieving high 11ft-to~drag (L/D) ratios at high transonic and 
low supersonic soeeds (Ref 1»2). Prominent among the features of this uncon- 
ventional aircraft is a wing of relatively large unswept aspect ratio (of the 
order of 8-12) which can be pivoted so that It presents Itself at an oblique 
angle to the flow (Fig. 1). Theoreti: «lly, this asymmetrical sweeping of a 
high aspect ratio wing results In t ' ry efficient wing shape in this speed 
range (Ref. 3). That this theory Is v:.11d has been shown in wind tunnel tests 
at the NASA Ames Research Center (Ref. 4). While demonstrably advantageous tc 
the aerodynamicist, such a design suggests potential structural stiffness and 
strength difficulties which deserve- careful consideration. Of these difficul- 
ties, the areas of static aeroelasticity and flutter behavior of the wing seem 
most important and worthy of analysis. 

The term static aeroelasticity is commonly applied to aeroelastic problems 
where inertia effects can be safely neglected. Control effectiveness and the 
redistribution of airloads on a flexible aircraft are prime examples of prob- 
Ic s where the equilibrium state of the flight vehicle Is highly dependent upon 
the interaction between the airloads and the flight structure. On the other 
hand, static wing divergence provides an example of a static aeroelastic stabi- 
lity problem. Since this latter problem is of no concern to a freely flying 
oblique winged aircraft (Ref. 2), the attention of this paper will be focused 
only on the problem of flexible airload redistribution and the lateral 
trim requirements of flexible wings which are asymmetrically swept. 

Earlier studies have focused attention of stability and control character- 
istics, which might be of importance to oblique wina aircraft design (Refs. 5, 
6,7). In Ref. 7, the author very briefly explored the possible Influence of 
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static aeroelasticity on oblique wing design. Further work In this area has 
) shown that the analysis presented In that reference can be extended and Inter- 
preted In such a way as to yield meaningful results. For this reason, the 
present study Is divided Into three parts: the analysis of the potential aile- 
ron control requirements for an oblique winged aircraft to ensure lateral 
equilibrium at various flight speeds; the presentation, analysis and comparison 
of an alternative mode of ensuring lateral equilibrium, the use of so-called 
built-in or geometric twist; and, finally, an assessment of the validity of the 
assumptions used and results obtained In the latter two studies. 

Discussion 

Some Insight into the magnitude and importance of the oblique wing lateral 
trim problem can be obtained through the study of the simplistic elastic wing 
model considered in Ref. 7. The analysis of this model, shown in Fig. 2, will 
seek to determine the extent to which aileron deflection, or some other means 
^ of control, is necessary to ensure lateral equilibrium. Although the analysis 
of the use of aileron deflection for this same model has been briefly discussed 
In Ref. 7, the essential features of that analysis will be reviewed here and 
the results extended because they bear heavily upon the ensuing geometric twist 
analysis. 

The idealization shown in Fig. 2 represents a continuous, uniform-property, 

elastic wing, of constant chord, swept at an angle A to the direction of 

flight. The wing is uncambered and has full- span ailerons for lateral control. 

For the formulation of the analytical model , the wing is assumed to behave as a 

slender beam with a straight elastic axis; torsic.ial flexibility of the wing is 

ignored to simplify the analysis. Finally, aerodynamic strip theory is used to 

describe both the applied loads^and the aeroelastic loads. With reference to 
\ * 

Fig. 2, the governing differential equation of equilibrium of the flexible 
wing, under the usual elementary beam theory assumptions can be written (in 
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nondimensional form (cf. Ref. 8, pp. 479-481) as: 


u 


d w , X dw 
— r + — 

dn dn 


El 


B6(n) 


-1 


< n < 1 


( 1 ) 


where w(n) = bending deformation, nondimensional i zed with respect to L. 

p^ = constant load per unit length in the n direction. 

El = bending stiffness, a constant. 

3 

X = qcc^.L sin A cos A/EI. 

B = 9CC|^gL^ cos^ A/EI. 

q = dynamic pressure. 

c, = two-d’:iensional sectional lift curve slope per unit chordwise 

angle of attack. 

c, . = two-dimensional sectional lift curve slope per unit aileron 
^ deflection. 

If this wing were to be clamped at its center, a roll moment would be gen- 
erated if ailerons were not applied in an anti symmetrical manner. This roll 
moment occurs because of the well-known tendency of sweptforward wings to 
develop lift forces when deflected upward and the tendency of sweptback wings 
to lose lift when similarly deflected upward. In the wind tunnel, a wing roll 
moment can be counteracted by the tunnel mount; in free flight, some mode of 
lateral control must be used. The analysis which follows assumes that longi- 
tudinal trim is supplied by control surfaces which do not enter this problem. 

As shown in Rer. 7, the full span ailerons may be used to null out the 
aeroelastic roll moment on the oblique wing if they are deflected such that the 
following relation is satisfied. 


o 




5(n) 


; 6 . 


- 6 , 


0 < n < 1 
-1 < 0^0 


( 2 ) 


\ 


Moment equilibrium about the roll axis requires that 



The assumption of the existence of a very stiff wing pivot support struc- 
ture at the aircraft centerline allows one to separate Eq. 1 into two regions, 
thus simplifying its solution. In this case, each wing half is assumed to be 
clamped at the aircraft centerline, or wing root position, and free of bending 
moments and shear at the wingtips. With these assumptions, a closed-form solu- 
tion for in Eq. 3 may be found. Because the term dw/dn enters into the 
aeroelastic load computation, this derivative is of more significance than the 
deflection w(n) itself. For this reason, the solution to Eq. 1 is usually pre- 
sented in terms of the variable r(n) where 




The variables and Tj^ are functions of the aeroelastic parameter X and are 
also to be found in the Appendix. 

It is found that, if the sweptforward wing portion is clamped at the root 
and in the absence of aileron application, the wing will undergo static diver- 
gence at a value of X= 6.33 (r=f. 8, 311-317). The value of the parameter x 
thus provides one measure of t'lc influence of static aeroelasticity in this 
type of problem. A relatec parameter which is sometimes useful is the variable 
q*. The parameter q* u .otes the ratio of the flight q to the divergence q, 

%ir 

q* = q/qQiv 

If the variables A, El, L and c are held fixed, then 

q*=X/6.33 (8) 

The magnitude of the aileron deflection parameter 6^ may b? conveniently 
examined if we look at the behavior of another parameter, y, found from a mani- 
pulation of Eq. 6. 

Y = 6.33 5^y' (p^L^/EI) ] /[ V ! 

A graph of y versus q* is given in Fig. 3. Two features of the relation shown 

in Fig. 3 are worthy of note. First of all, as q* approaches the value 4.335, 

the parameter y tends to an infinite value because, at this value of q*, T^^ = 

T„. Also, it is seen that the value of y is nearly equal to unity for values 
K \ 

of q* in the range 0 < q* < 1.50. If this latter observation is used, together 
with Eq. 9, then 

\ 
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(10a) 


or 6^ (E^’) te jtan A ^ 6.33 (10b) 

Eq. 10b can be rendered more meaningful if is related to the Ig flight 
condition of an aircraft of gross weight W. represents the load per unit 
length (assumed constant) along the swept span caused by the airloads on the 
rigid wing; a relation between p^ and W is readily obtained if we ignore both 
the aileron induced airloads and the aeroelastically induced airloads. This 
relation is 

2P„L 1 W (11) 

Therefore, the first coefficient in Eq. 10b becomes 

n, WL.2 (12) 

El " 2EI 

The degree to which Eq. 12 is an approximation will be discussed later. How- 
ever, it should be remarked here that the inclusion of aeroelastic lift into 
the analysis causes the term on the left hand side of Eq. 12 to be a fraction 
of the term on the right. This fraction is very close to unity for reasonable 
values of q*. By substitution of Eq. 12 into Eq. 10b, one obtains the result: 

The most striking feature of Eq. 13 is that 6 ^ does not depend on the dynamic 
pressure and therefore is not an explicit function of the flight speed. The 
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WL 

El 


n r 

J L 


-Lg tan A 
'L6 


H2.66 (radians) 


(13) 


- 6 - 


parameter does, however, depend upon the wing flexibility and the sweep 
angle A. Although formulated for a Ig flight condition, the weight W enters 
the equation linearly and could just as well have been written as MW for an Ng 
condition. This feature of Eq. 13 means that, once the aircraft Is trlmned for 
one flight speed at a given sweep angle A, It Is trimmed for all flight speeds 
at that same sweep angle and altitude. To Illustrate the order of magnitude of 
the aileron deflection which might occur for a large transport aircraft, let us 
use the following data which are chosen to be typical of this class of aircraft. 

W = 400,000 lb. 

El = 20.0 X 10^^ lb-in. 

^L6 

L = 1200 in. 

For this data, we obtain, from Eq. 13: 

= 3.26 tan A J (degrees) (14) 

Mt 45“ of sweep, a full-span asymmetrical aileron deflection of 3.26“ is neces- 
sary to ensure lateral equilibrium. Although not extremely large, such a 
deflection might have an adverse effect on aircraft yaw trim and cruise L/D. 

Turning to the second objective of this paper, the previous discussion has 
illustrated but one method of controlling the tendency of the oblique wing to 
develop a rolling moment in flight. One method of improving aerodynamic per- 
formance in aircraft has been the use of geometric angle of attack distributions 
or "built-in twist," For this application, the flight structure Is geometric- 
ally tailored to satisfy some performance objective. 

Wind tunnel tests of "rigid" oblique wings and the analysis of such wings 
by methods employing accurate aerodynamic theories have shown that there Is a 
tendency for the lift distribution to build up on the sweptback wing. This 


tendency causes a roll moment opoosite in direction to that caused by aeroelas- 
tic effects. To cancel out this adverse situation, it has been proposed 
that some amount of upward wing geometric curvature or dihedral be used to 
alleviate this roll moment (Ref. 9). 

To understand how a swept wing with a built-in deflection can develop an 
angle of attack with respect to the airstream, consider the situation shown in 
Fig. 4 (the reader is referred to a similar, more complete discussion given in 
Ref. 10, pp. 474-479). Since small rotations transform as vector quantities, 
it is seen that the angle of attack «, seen by the freestream parallel to the 
flight direction is, for a sweptback wing section, given by the expression 

« = 0 cos A - t|) sin A (15) 

For a sweptforward wing, the above relation is modified by substitution of a 
"plus" sign for the negative sign before the second term on the right hand side. 

If 0, the twist angle along the swept axis, is zero then the inclination 
of a sweptback wing such that the swept axis forms an angle ij; with the horizon- 
tal plane, i.e. a dihedral, results in a constant negative angle of attack 
along the wing. Conversely, a dihedral on a sweptforward wing results in a 
constant positive angle of attack. Thus, for an oblique wing, a built-in dihe- 
dral angle tends to generate a set of lift forces which roll the aircraft so as 
to elevate the sweptforward wing tip; an anhedral would produce just the oppo- 
sit.^ effect. The key point here is that the anhedral -dihedral effect is, in 
esr>r>nce, a built-in geometric twist effect. The purpose of the ensuing discus- 
sion is to show how this effect can provide lateral equilibrium without any 
aileron action. 

^ith the assumptions used jn the previous aileron study and once again 
using chorduise Cross sections, the nondimensional governing equation of equi- 
librium for the flexible wing with an initial built-in dihedral angle 

\ 

- 8 - 






distribution ii(n) (i»(n) is the angle formed by the swept axis and the hori 
zontal plane) is as follov/s: 


d^w ^ ^ _ ^0^ - X4»(n) i ^ 06) 

dn El 


The definitions of the quantities other than i(i(n) used in Eq. 16 are identical 
to those used in Eq. 1. To cancel the aeroelastic roll moment, a constant ini- 
tial built-in anhedral angle is used such that i|»(n) is defined as: 


'(»(n) 


+4; 

^0 

-1 < n 0 

o 

1 

O 

A 

1 A 


(17) 


The comparison of Eqs. 16 and 17 with Eqs. 1 and 2 shows that they are 
made identical if 


=6 

r 




(18) 


To guarantee roll moment equilibrium, must satisfy the relation 


= f n^dn*/* n r(n) dn 
J -1 dn J -1 


(19) 


The solution for the variable r(n) = ^ in this problem is identical to that 

dn 


presented for the previous aileron problem if the term replaces the term 


in the expressions for T|^ and given in the Appendix. Substitution of 


these expressions into Eq. 19 and subsequent integration yields a relation for 


given by: 
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X 

El 

— mm 




( 20 ) 






The expressions for Tj^ and Tj^ are those presented in the Appendix. 

Assumptions similar to those given for the aileron analysis lead us to an 
approximation for Eq. 20 given by: 

♦o'— (2’> 

° [rr J 

An examination of Eq. 21 and comparison of this relation with Eq. 13 
reveals how remarkably efficient the use of built-in twist, in the form of the 
initial anhedral, is for the oblique wing roll equilibrium problem. Eq. 21 has 
no factor tan A, reflecting the fact that the same airfoil sections 

which are causing the roll equilibrium problem are also being used for its 
solution. As the influence of aeroelasticity increases with the angle of sweep, 
so too does the counteracting built-in twist effect, since the sine of A in Eq. 15 
increases with A. 

From Eq. 21, it is seen that the anhedral angle of the swept axis is not a 
function of flight speed. Use of the same typical parameters as were consid- 
ered for the aileron example results in a value for of 

= 1.30 degrees (22) 

This initial anhedral angle is small and corresponds to a situation where the 
wing tips are initially located a distance of 27.2 inches below a horizontal 
plane passing through the fuselage centerline, a rather small distance when 
compared to the 100 foot semi -span. 

Since the simplistic formulas given by Eqs. 13 and 21 are approximations, 
even for the idealized model used in this analysis, two additional tasks remain 
before this study can be considered complete. The first of these tasks is the 
detei^ination of the extent to Which the aeroelastic lift contribution affects 
the accuracy of these strip theory expressions. The second task is to 
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determine to what extent the use of aerodynamic strip theory to describe the 
loads affects the problem solution. For the sake of brevity, we will examine 
only the geometric twist problem, since the solutions to this problem and to 
the aileron problem are mathematically similar. 

To accomplish the first task, we begin by summing all of the vertical 
forces which act upon the idealized wing and the equate these forces to the 
aircraft gross weight; the following expression results. 


2Po'-^ 

rr 




dw dn _ WL^ 

^ - IT 


(23) 


The integral term in Eq. 23 represents the relative vertical distance between 
the wing tips. Physical reasoning leads one to expect that this will be a 
small negative quantity. It is this latter term which is ignored in the formu- 
lation of Eqs. 12 and 20. 

Direct integration of the integral term in Eq. 23 using the expressions 
for r(n) and subsequent solution for yields the relation: 


El" 


WL^ 

2EI 




(24) 


U, and U„ are defined in the Appendix. The substitution of Eq. 24 into Eq. 20 

L K ^ 

gives the exact .solution for 


(25) 


= 1 



WL^ 

2X 

_ 


IT 


The coefficient of the factor j^WL^/Elj is a function of the aeroelastic 
parameter x. Eq. 21 approximates the value of this coefficient to be 1/12.66 
or 0.07900. This approximation is the result of taking the value of y to be 


unity and ignoring th'' aeroelastic lift. The ratio ^^/(WL /El) is shown in 
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precluded. Even more Important 1s the realization that the use of strip theory 
results in the overestimation of loads towards the wing tips. 

To investigate these possible shortcomings an analysis, based upon the 
theoretical model developed in Ref. 11, was carried out. The aerodynamic theory 
used in Ref. 11 is based upon a modification of the Weissinger-L method 
detailed in Ref. 12. The structural model assumes the wing to be a beam with 
bending-torsion flexibility. The analysis method detailed in Ref. 12 Is a 
matrix method which permits spanwise variations in wing elastic properties; its 
chief advantages are the ease with which it may be programmed for the computer 
and the relatively short computer run times needed for static aeroelastic 
analysis. 

To test the validity of the simplistic relations derived previously in 
this paper, a computer study using the model just described was conducted on an 
elastic wing model similar in size to a set of elastic wings used in wind tun- 
nel tests at Virginia Polytechnic Institute (Ref. 13). The idealized model has 
the following properties: 
c =4 inches 
L =20 inches 
c, =6.28 per radian 
El = 1000.0 lb-in. 

The torsional stiffness input to the computer program was chosen to be ten 
times the bending stiffness so that only bending flexibility was important to 
the analysis. The results of the analysis of this idealized, uncambered wing 
give the value of the initial built-in anhedral in terms of the wing loading 
parameter, W/S Ib/ft . Given the scale of the idealized wing model, a value of 
W/S of the order of 1 to 2 psf.tis probably representative when compared with 
the large transport data given previously. For this small model, the use of 
Eq. 21 givfi5 an anhedral value of 

-13- 


degrees 


(27) 


^ 2.012 X. 



The above mentioned computer analysis was run at a Mach number of zero and 

■ / ’>■ 
at three sweep angles, A = 15®, 30® and 45®. The studies were run for q* 

' : : ‘ I : 

values in the range 0 < q < 2.5. Attention is called to the fact that q^jy 
changes with A. The results of this study are displayed in Fig. 6. 

In Fig. 6, the required anhedral angle (per unit W/S psf.) is plotted ver- 
sus q* for the three sweep angles considei'ed. A solid straight line is drawn 
to represent the strip theory prediction. From this figure., it is seen that at 
each va:iue of A, the required value of is a function of q*. The change in 
required is dramatic in the range 0 < q* < 0.50, but less so above q* = 0.50. 
In all cases the strip theory formula overestimates the required value of 
the reason for this overestimation is that the use of strip theory results in 
an ideali-Zation that is too flexible from a static aeroelastic standpoint. 

That this latter observation is true can be seen in the fact that strip theory 
underestimates the value of q^jy by about 30% when compared to the numerical 
method being used here. 

Of further interest is the observation that, at low values of q*, a nega- 
tive value of tjj is required. This corresponds to the observation 

in Ref. 9 that some slight upward curvature of the wing is necessary to correct 
a small roll moment caused by the lift distribution shifting toward the down- 
stream wing tip. These small values of q* correspond to what the author would 
term a "rigid" wing. 

Some potential importance may be attached to the fact that crossover 
points exist in the curves presented in Fig. 6 where the required anhedral 
changes sign. It is to be noted that these changes in sign occur at progfes- ^ 

. \ ^ ^ i 

sively larger values of q* as A increases. At this point, by definition, no 
lateral trim is required to maintain equilibrium. This occurs because the roll 



moment caused by the airload buildup towards the sweptback wing tip Is just 
cancelled by the aeroelastica‘'iy induced load buildup towards the sweptforwJird 
wing tip. On an actual aircraft, this crossover or zero point would be a func- 
tion of such parameters as wing planform shape, sweep angle, Mach number and 
stiffness distribution. 

At the crossover point, no built-in anhedral or aileron action is neces- 
sary, no matter what the wing loading or load factor. The desirability of 
designing the wing so that the zero point occurs as closely to the cruise speed 
as possible, if not precisely at this speed, is obvious. Since so many other 
design objectives must be met by the structural engineer, this latter objective 
may be difficult to fulfill in practice. 

Conclusions 

Some unique static aeroelastic problems posed by the asymmetrical sweeping 
of a high aspect ratio wing have been examined through the use of a simplified 
model. Although the use of such an idealization to represent problems which 
are likely to occur in actual aircraft designs and the solutions to these 
problems is an oversimplification, the author nevertheless believes that 
several useful conclusions or guidelines may be drawn from this 'study. 

The use of strip theory to represent aerodynamic and aeroelastic loads 

leads to answers which overestimate the amount of aileron input or geometric 

twist necessary to ensure lateral equilibrium, particularly at low values of 
★ 

q . In addition, the error introduced by strip theory increases as A in- 
creases. The simplistic formulas derived from strip theory assumptions thus 
give\:onservativ8 results muchas they do when' used in conventional static 
aeroelastic stability studies. 


The results derived through the use of irore accurate aerodynamic theory 
show that as the parameter q* increases, increased demand is made on the method 
which guarantees lateral equilibrium. For small values of q* near zero, few 
demands are made to guarantee lateral equilibrium. The parameter q* can be 
made small, given a constant operating q, by increasing the clamped divergence 
q of the aircraft by one of several conventional methods; these methods include 
stiffening the structure or redistributing the wing area so that more of the 
area is inboard. Fo»^ aircraft which are designed in a conventional manner to 
the usual strength and stiffness criteria, the amount of modification to pre- 
clude the roll problem discussed in this paper is probably minimal. The normal 
"droop" of an aircraft wing due to gravity provides some anhedral effect. 

Of potential theoretical interest in the areas discussed in this paper is 
the use of structural modifications to further improve static aeroelastic 
performance. Modifications such as asymmetrical wing stiffening or redistri- 
bution of wing stiffness to bring the crossover or zero point near to the cruise 
q might prove to be worti'ishile. A similar study of the use of various wing 
planforms and their relative merits might also be in order. 

To summarize, this problem of asymmetrical wing static aeroelastic equili- 
brium is one which must certainly be considered by the designers of such an 
aircraft. It is likely, however, that after all the conventional design cri- 
teria are met, this additional unique criteria will cause few, if any, 
additional problems for this aircraft configuration. 
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TABLE I 



q* 

,|;q/(WL^/EI) 

(Exact value) 

Per cent error in 
approximate value 

(PqL^/EI)/(WL^/ 2EI) 
(Exact value) 


0.00 

0.0750 

+ 5.33 

1.00 

0.25 

0.0751 

+ 5.19 

0.998 

0.50 

0.0753 

4.91 

0.991 

0.75 

0.0756 

+ 4.50 

0.980 

1.00 

0.0761 

+ 3.81 

0.964 

1.50 

0.0776 

+ 1.81 

0.916 

2.00 

0.0797 

- 0.878 

0.846 

2.50 

0.0827 

- 4.47 

0.749 

3.00 

0.0867 

- 8.88 

0.618 

3.50 

0.0921 

-14.2 

0.444 

4.00 

0.0994 

-20.5 

0.207 ' 


The approximate value of t|)^/(WL^/EI) is 0.0790. 



APPENDIX 


The analytical expressions for rj^(n)» rj^{n), T^(a), Tj^(a), Uj^(x) and 
Uf,(x) used in the body of this paper are presented below. In the region 
-1 < n i 0, r(n) is given by: 


r,(n) = ^ - M„] * 

, g-ad+n) + 2e^(l+n)/2 

I _ 


..-a i f 

e + 2e cui t 


where a 


J/3 


and f = a(3^''^/2). 


In the region 0 < n l 1, r(n) is given by: 
3 


PrU) = 


i_r?oi£ 

aHEI 


+ p6. 


. e^O-n) + f(l-n) 

1 “ ; rTTo 


e® + 2e"®'^^cos f 


The expressions for Tj^ and Tj^ are given by: 


e’3a/2 _ cos f + 3^^^sin f 






a^(e 

-3a/2 . 

*■ 2cos f) 

e3a/2 , 

cos f 

- 3^/2 f 

a^(e 

3a/ 2 ^ 

2cos f) 

U|^ and 

U|^ are 

; 

cos f + 

J/2 • f -3a/2 

3 ' sin f - e 

a(e 

-3a/2 

+ 2cos f) 

cos f - 

,1/2 . ^ 3a/2 

3 'sin f - e 


\ 


^ + 2cos f) 


(Al) 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


\ 
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